COMBINATORIAL MODELS OF RIGIDITY AND 
RENORMALIZATION 



J. BARRE 

Abstract. We first introduce the percolation problems associ- 
ated with the graph theoretical concepts of (k, Z)-sparsity, and 
make contact with the physical concepts of ordinary and rigid- 
ity percolation. We then devise a renormalization transformation 
for (fc, Z)-percolation problems, and investigate its domain of va- 
lidity. In particular, we show that it allows an exact solution of 
(k, ^-percolation problems on hierarchical graphs, for k < I < 2k. 
We introduce and solve by renormalization such a model, which 
has the interesting feature of showing both ordinary percolation 
and rigidity percolation phase transitions, depending on the values 
of the parameters. 



1. Introduction 

Consider an ensemble of bars glued together at joints, around which 
they can freely rotate, embedded in a two dimensional space. The 
ensemble of bars and joints is naturally associated with a graph: joints 
are the vertices of the graph, and bars its edges. The question of 
whether or not this ensemble constitutes a rigid body is a fascinating 
one. It found a convenient answer with Laman's theorem pQ, which 
characterizes rigidity in a purely graph theoretical way[] It implies for 
instance that a graph with N vertices must have at least M = 2N — 3 
edges to be rigid. Here the number 2 corresponds to the number of 
degrees of freedom of a joint and the number 3 to the number of degrees 
of freedom of a rigid body, in two dimensions. The concepts of (k, /)- 
sparsity and (k, /)-tightness, introduced a long time ago in graph theory 
[2], allow to generalize the problem [3]. For instance, bar-joint rigidity 
in two dimensions corresponds to (k, I) = (2,3). Some values of (k, I) 
have a clear physical meaning in terms of other models of rigidity, 
others do not, see Sec. [2] 



Throughout this paper, we will be concerned only with generic rigidity, which 
can be reduced to a combinatorial problem. See for instance [5T] for a discussion 
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The percolation problem associated with bar-joint rigidity has drawn 
the attention of physicist's for a long time jl]. It can be stated as fol- 
lows: take a large lattice, where each site is a joint and each bond 
a bar; consider that each bar is effectively present with probability p. 
When p is increased, the system goes from a "floppy" phase where only 
small rigid clusters exist, to a "rigid" phase, where there is one macro- 
scopic rigid cluster. In between lies the rigidity percolation transition. 
To study this problem, Jacobs, Thorpe, Duxbury and Moukarzel have 
introduced the so called "pebble game" algorithm, which implements 
in an efficient way the combinatorial characterization of Laman [HI E], 
see also [9] . The interpretation of the numerical simulations on regular 
two dimensionnal lattices have fueled some debate, but the following 
picture seems favored: the rigidity percolation transition in such 2D 
bar-joints lattices is second order, in a universality class different from 
ordinary connectivity percolation [HI US] • We note that rigidity perco- 
lation may also be studied on trees or random graphs, in which case it 
is exactly soluble, and displays a first order transition, at variance with 
ordinary percolation [TT]. To our knowledge, there exists no analytical 
theory, even approximate, for the critical exponents of 2D bar-joint 
rigidity percolation on regular lattices. 

Except for some special cases, such as ordinary percolation and the 
bar-joint rigidity percolation described above, percolation problems as- 
sociated to (k, I) sparsity seem to have received relatively little atten- 
tion in the physics literature (see Sec J2] for a more detailed discussion). 
In particular, little is known about the associated percolation transi- 
tions and their universality classes. A natural idea to get insight into 
this problem is renormalization theory. The purpose of this paper is 
to introduce a renormalization tool to study (k, Z)-percolation, gener- 
alizing the method of [16] . In Sec. |3j we will determine the conditions 
under which the renormalization procedure is useful. In particular, 
we will show that it allows a complete solution for some (k, /)-rigidity 
models on hierarchical graphs [HI [16]. In Secjlj we present and solve 
explicitly a model of (3, 3) percolation (which has a physical interpre- 
tation as bodies and bars percolation in two dimensions). We will show 
that the renormalization transformation on this model has two trivial 
fixed points, corresponding to the full and empty graph respectively, 
and two critical fixed points. One fixed point, corresponding to ordi- 
nary percolation has three unstable direction, whereas the other one 
has one unstable direction, and thus governs the critical behavior for 
almost all values of the parameters. This model provides an example of 
a system showing both ordinary and rigidity percolation transistions. 
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2. Rigidity and graph theory 

We first briefly recall the graph theoretic description of bar-joints 
rigidity [I]. The basic idea is constraint counting: each joint has two 
translational degrees of freedom, and each bar, by fixing the distance 
between two joints, removes one degree of freedom. Since a rigid body 
in two dimensions has three degrees of freedom (two translations and 
one rotation), an ensemble with N joints needs at least M = 2N — 3 
bonds to be rigid. In addition, these bars must be "well distributed" 
among the joints: if there are exactly 2N — 3 bonds, there should be no 
subensemble with n joints containing strictly more than 2n — 3 bonds. 

Clearly, the numbers 2 and 3 in this description are adapted to bar- 
joints rigidity theory in two dimensions, but it is natural to generalize 
the idea to other pairs of numbers (k,l). [20| [3]. It is also natural to 
allow for multiple edges in a graph: such graphs are called multigraphs. 
Although, to our knowledge, they have not appeared in this form in the 
physics literature, the following definitions about (k, Z)-rigidity are not 
new [3J; we recall them for self-consistency of the article. Throughout 
the article, we shall always consider multigraphs, although we shall 
keep calling them graphs. 

Given a graph G = (V, E), where V represents the set of vertices of 
G and E the set of its edges, we say that G' = (V, E') is a subgraph 
of G if V C V, E' C E, and the edges in E' connect vertices in V . 

Definition 2.1. For k > 1 and < I < 2k — 1, a multigraph G with N 
vertices and M edges is said to be 

- (k, /)-sparse if every subgraph G' of G with n > 1 vertices contains 
at most m = kn — I edges; 

- (k, /)-tight if in addition M = kN - I; 

- (k, Z)-spanning if there is a (k, Z)-tight subgraph spanning all vertices 
ofG; 

- (k, Z)-redundant if it is not (k, /)-sparse. 

Two vertices A and B are said to be (k, /)-rigidly connected if there is 
a (k, /)-tight subgraph containing A and B. 

When there is no ambiguity, we shall sometimes omit the (k, I) in- 
dices, and talk about sparse, tight, and redundant graphs. 

Notice that if I > 2k, these definitions make no sense: there is no non 
trivial (k, Z)-sparse graph in this case, since the simplest graph with 2 
vertices and 1 edge is not (k, Z)-sparse. More generally, the edges in a 
(k, /)-sparse multigraph have multiplicity at most 2k — I. 

These graph-theoretic definitions may be immediately translated in 
the rigidity vocabulary. A (k, Z)-tight graph is a (k, I) isostatically rigid, 
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or minimally rigid, graph: it has just the right number of constraints 
to be rigid, and removing any constraint destroys rigidity. A (k, l)- 
spanning graph is a rigid graph according to the (k,l) rule. 

Redundant constraints play an important role in bar and joint rigid- 
ity theory. It is straightforward to generalize the concept to (k, l)- 
rigidity. A graph possesses one or more redundant edge if and only if 
it is not sparse, ie if it contains a subgraph with n vertices, m edges, 
and m > kn — I. To define and count the number of (k, Z)-redundant 
edges in a graph G = (V, E), one may proceed as follows: 

(1) Start with the empty graph (V,0), without edges 

(2) Add the edges of E one at a time 

(3) Check if this addition creates a redundant subgraph, with n 
vertices, m edges, and m = kn — I + 1. 

(4) If this is the case, then reject the newly added edge and add 
one to the count of redundant constraints 

(5) Continue until all edges in E have been added. 

This procedure actually follows the steps of the algorithms used to 
study (2, 3)-rigidity 0, EE 13 [10] , and more generally (k, /)-rigidity [3]. 
A priori, the number of redundant edges computed in this way could 
depend on the order chosen to add the edges. It is not the case, so 
that the number of redundant edges is a well-defined concept. This is 
also the minimum number of edges that one has to withdraw to make 
a redundant graph sparse. In general, there is a freedom of choice in 
the edges to be removed to turn a redundant graph into a sparse one, 
but their number is fixed. 

Beyond bar-joint 2D rigidity, which is a physical realization of (2, 3)- 
rigidity, (k, I) -rigidity encompasses several known situations: 

• (1, l)-rigidity corresponds to ordinary percolation, in any di- 
mension. 

• (3, 3)-rigidity corresponds to the rigidity of systems composed 
of rigid bodies and bars in 2 dimensions. 

• More generally, (fc, fc)-rigidity is related to rigidity of bodies and 
bars systems in higher dimension [20J . 

• (3, 6)-rigidity is related to the rigidity of bond-bending networks 
in 3 dimensions [2T| 122] . In this case 1 = 6 = 2k, which is outside 
the usual range; the definition of sparsity needs to be slightly 
adapted. 

• (2, 0)-rigidity corresponds to the graph theoretic concept of 2- 
orientability, which has been used to study bar-joint rigidity of 
random networks |23j . 
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Clearly, a percolation problem may be associated to each case, by mon- 
itoring the largest (k, /)-rigid cluster in a graph. Understanding the 
general properties of these percolation processes, such as the order of 
the transition, the universality class, etc... is a largely open question, 
although several cases have been studied in details: among those, or- 
dinary percolation is of course the best known by far. Percolation for 
bond-bending networks has been studied, because of its relevance in 
modeling network glasses and proteins [12], but we are not aware of 
studies on critical exponents. In [13], C. Moukarzel studies percolation 
problems on trees and random graphs for various numbers k (g in the 
notations of [13])[^J It is also known that when some bars and joints 
form a rigid component, they behave in every aspect like a rigid body. 
This idea has actually been numerically exploited [14J. Then one may 
expect that body-bar ((k,l) = (3,3)) and bar-joint ((k,l) = (2,3)) 
rigidity percolation belong to the same universality class [To] . 

Summarizing these informations, one may identify several conjec- 
tures and questions: 

(1) (1, 1) percolation (that is ordinary percolation) and (2, 3) perco- 
lation on 2D regular lattices belong to two different universality 
classes: this conjecture seems well documented numerically. 

(2) One may also conjecture that (2, 3)-percolation (bar-joint) and 
(3, 3)-percolation (body-bar) on 2D regular lattices belong to 
the same universality class. 

(3) On random graphs, it is known that the threshold for (2, 3)- 
percolation is identical to the threshold for (2, 0)-percolation [23] 
One may then wonder if, on random graphs, the threshold for 
(k, /)-percolation could be independent of /, for < I < 2k — 1. 

A more general question is: given for instance a regular lattice in d 
dimensions; what does determine the order of the transition and the 
universality class of (k, I) percolation? Efficient algorithms, generaliz- 
ing to (k, Z)-sparsity the original pebble-game devoted to (2, 3) sparsity, 
have been recently introduced [3J, so that a numerical investigation of 
this question is in principle possible. To our knowledge, it has not 
been undertaken. A natural analytical tool to progress in the studies 
of these percolation processes is the renormalization group. In section 
[3] we introduce a renormalization transformation, and we prove that in 
some sense it is equivalent to study the rigidity properties of the orig- 
inal and renormalized graphs, when the parameter range is restricted 



As a tree alone is never rigid, the leafs are usually glued to a busbar to study 
rigidity for trees. 
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to k < I. This allows to solve exactly models on hierarchical graphs for 
k < I; in Sec. [4], we solve such a model for k = £ = 3. 

3. Renormalization rule for a graph 

The renormalization rule we are about to describe is a generalization 
of the rule introduced on [T5] for (k, I) = (2, 3). In this case, this rule is 
intuitive enough so that it does not really require a justification relying 
on graph theory, see [IE] (this rule has been rediscovered independently 
in [Ej). Things are not so obvious for general (k,l), and in particular 
we have to impose the restriction k < I. In order to be sure to make 
correct statements, we feel it is useful to give precise graph theoretical 
definitions, and to proceed rigorously, not relying only on intuition. 
This certainly makes the article heavier to read, but is necessary to 
identify for instance the restriction k < I, which is not completely 
intuitive. 

From now on, we assume k < I < 2k. We will denote \F\ the cardinal 
of an ensemble F. 

3.1. Renormalization of a whole graph. Consider a multigraph 
H = (V, E), where V is the set of vertices and E the set of edges. We 
singularize two sites A,B e V. We would like to replace the graph H 
by an "equivalent" (in a certain sense) graph with only two vertices 
A and B, and a certain number of edges Mab linking A and B. We 
define this renormalization step as follows: 

Let us call Sab the ensemble of all (k, /)-sparse subgraphs of H con- 
taining A and B. For h = (Vh,Eh) in Sab, the number of its vertices 
and edges is respectively \Vh\ and \Eh\. We define 

r AB {H) = max [\E h \ - k\V h \] 
heSAB 

Note that since h is (k, /)-sparse, tab^H) < —I. Then Mas is given by 

(3.1) Mab — max (r ab(H) + 2k, 0) 

We have < Mab — 2k — I. The renormalized graph, which contains 
vertices A and B and Mab edges, is sparse; it is tight if and only if A 
and B are (k, /)-rigidly connected in the original graph H. To better 



understand formula (3.1), let us explain in details what it means for 
some particular (k,l). 

Ordinary percolation 
(k, I) = (1, 1) corresponds to the ordinary percolation case. Then 2k — 
1 = 1, and the number of renormalized edges between A and B is or 
1. Sparse graphs in this case are graphs without any loop. One sees 
that a loopless subgraph h = (Vh,E h ) with \E h \ > \ V/J — 1 is necessary 
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Figure 1. Example of renormalization of a simple 
graph with 4 vertices, k = I = 1 (ordinary percolation). 
The upper graph is renormalized into one edge; the lower 
graph is renormalized into the empty graph. Note that 
both graphs contain one redundant edge. 




FIGURE 2. Example of renormalization of a simple 
graph with 4 vertices, k = 2 and / = 3 (bar-joint rigidity 
percolation). The upper graph is renormalized into one 
edge, as A and B are rigidly connected; the lower graph 
is renormalized into the empty graph, as A and B are 
not rigidly connected. 



for "percolation" across the graph. On the elementary "Wheatstone 
bridge" graph of Fig. |3.1[ it is easy to see that this coincides with the 
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FIGURE 3. Example of renormalization of a simple 
graph with 4 vertices, k = 3 and I = 3 (body-bar rigidity 
percolation). The upper graph is renormalized into three 
edg GS, clS the bodies A and B are rigidly connected; the 
lower graph is renormalized into two edges only. Note 
that the upper graph contains one redundant edge, for 
instance the dashed one. 



usual renormalization prescription: replace the graph by an edge if and 
only if A and B are connected. 

Bar-joint rigidity percolation 
In this case (k, I) = (2, 3). Again 2k — I — 1, and the number of renor- 
malized edges between A and B is or 1. The rule says in this case: 
"Renormalize the graph by a edge if and only if A and B are rigidly 
connected (ie belong to the same rigid cluster)". This is intuitive, and 
this is the rule which has been used in [IB]. Examples are given on 



Fig. 3.1 



Body-bar rigidity percolation 
In this case (k, I) = (3, 3). 2k — I = 3, so that the number of renormal- 
ized edges between A and B is 0, 1, 2 or 3. A and B have to be seen as 
rigid bodies, with three degrees of freedom each, and three edges are 
necessary for a rigid connection between them. Examples are given on 



Fig. 3.1 



The following property is useful to effectively compute renormalized 
graphs. 

Proposition 3.1. Let H = (V, E) be sparse, with A, B e V . Let e be 

an edge such that H' = (V,EU {e}) is not sparse. Then rAB^H) = 
r AB {H'). 



As a consequence of this proposition, the renormalizations of H and 
H' contain the same number of edges. In other words, adding redun- 
dant edges to a graph H does not modify its renormalization. We will 
use this as a tool to compute renormalized graphs in Sec. |4j 

Proof: Clearly, v A b(H) < tab{H'), since an element of Sab{H) is 
also in a natural way an element of Sab{H'). 

We call u and v the two vertices connected by the edge e; u, v G V. 

Let us now consider g' = (V',E') an element of Sab{H'), such that 
the maximum is attained in the definition of tab(H'). Thus \E'\ — 
k\V'\ = tab(H'). We want to construct a subgraph g" = (V",E") in 
S A b{H) such that \E"\ — k\V"\ > r AB {H'). This will prove r AB {H) > 
tab{H'). We may assume that e G E', otherwise it is enough to take 
g" = g' . Thus the vertices u and v are in V. 

Let G' = (Vq, Eq U {e}) be a minimal redundant subgraph of H' . V 
contains u and v the vertices connected by e, and G = (Vq, Eq) is a 
tight subgraph of H. We define V" = V U V and g" = (V",E") as 
the subgraph of H induced by the vertices V U Vo- g" is an element of 
Sab{H). g" does not contain e, but we will show that \E"\ — k\V"\ > 
r{H'). 

Assume first Vo C V. In this case V" = V. Then there is one edge 
eo in Eq U {e} that is not in E', otherwise g' would not be sparse, eo 
is different from e, thus eo is also in E". We conclude \E"\ — k\V"\ > 
\E'\-k\V'\=r A B{H'). 

We assume now that V \V' is not empty. We have the relations: 

V" = V'U (V \V) , E" D (E'\{e}) U (E \E') 
where U is a disjoint union. Indeed, e G E' and e ^ Thus 

(3.2) IV"! = \V'\ + \V Q \V')\ , \E"\ > \E'\ - 1 + \E \E'\ 

Let Gi = (Vq fl V, -Eo H E') the subgraph defined by the intersection 
of V and Vq. G\ is a subgraph of H, so it is sparse; it is also a 
proper subgraph of Gq, and it contains vertices u and v, so it is not 
tight, otherwise it would contradict the minimality of G' : indeed, (VbH 
V, E fl £" U {e}) would be a redundant subgraph of if', smaller than 
G' . We conclude that 

(3.3) \e n £'| < k\v nv'\-l 

We use now V = (V n V) U (Vb\V') and E = {Eq n £") U (Eo\E?). 
Since G is tight, we have 

(3.4) |£ n + \E Q \E'\ - k\V n V'\ - k\V \V'\ -1 = 



Using (3.3), we get 

\E \E'\ - k\V \V'\ > 
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Figure 4. On the left is schematized a graph G; the 
subgraph H satisfies the renormalizability condition. On 
the right, H has been renormalized into one edge, the rest 
of the graph is not modified. 



Putting this together with (3.2), we have finally 

\E"\-k\V"\ > \E , \-k\V'\-l+\Eo\E , \-k\V Q \V , \ > \E'\-k\V'\ = r AB (H') 

□ 

3.2. Renormalization of a subgraph. We have seen how to replace 
a graph by a certain number of renormalized edges. Now, we will see 
why and when it is licit to perform this renormalization step on a 
subgraph of a bigger graph. Consider a graph G = (V, E), Vh C V 
and H = (Vh, Eh) the subgraph of G induced by the vertices Vh- We 
will say that H satisfies the renormalizability condition if there exist 
two vertices A and B such that all edges in E linking H to the rest of 



the big graph G are connected to A or B; see Fig, 3.2 In the following, 
we always assume this condition is satisfied. 

For any subgraph of g C G, we can define the renormalized graph 
R(g), as follows: the part of g which is a subgraph of H is renormalized 



according to the rule detailed in 3.1; the remaining part of g is left 



unchanged. Note that if A or B, or both, do not belong to g, then the 
renormalized part of g contains no edge. 

Let us call G = R(G) = (V, E) the graph G where the subgraph H 
has been renormalized. We would like to show that it is in some sense 
equivalent to study (k, /)-percolation on G and G. For this purpose, 
we will prove the following: 

Proposition 3.2. In the setting described above: 

i) A and B are (k, I) -rigidly connected in G if and only if they are 
(k , I) -rigidly connected in G. 

ii) The number of redundant edges in G and G are related through the 
formula 

n re d(G) = n red (G) + n red (H) 



11 



where we have written n re d(F) for the number of redundant edges in a 
graph F. 

This last formula allows then to replace the problem of computing 
the number of redundant edges in G by the the problem of computing 
the number of redundant edges in the renormalized graph R(G) = G. 



Before proving 3.2 i) and ii), we state and prove a few lemmas. 



Lemma 3.3. Let G = (V, E) and R(G) = (V, E) be as above. Let 
g = (v,e) be a subgraph of R(G), and let us call h its renormalized 
part; h contains fh edges. Then there exists g subgraph of G, such that 

i) ~g = R{g) 

ii) If A or B do not belong to v, then, calling h = (Vh,Eh) the inter- 
section of g with H, E h is empty (ie h has no edge). 

Hi) If A and B belong to v, then h = (Vh,Eh) is sparse and satisfies 

fh = \E h \ - k\V h \ + 2k 

Proof: Note first that h has at most two vertices, A and B. If A or 
B do not belong to g, or if fh = then we choose for the vertices of h 
the same vertices as h, and put no edge. 

If fh > 0, we take h = (Vh,Eh) an element of Sab(H) such that 
\Eh\ — k\Vh\ + 2k = fh, and h minimal in the sense that no proper 
subgraph of h has the same properties. Then the renormalization of 
h produces exactly fh edges between A and B. Indeed, suppose it 
produces m' > fh such edges; then there exists h' = (V^, E' h ) subgraph 
of h, in Sab(H), such that \E' h \ — fc|V^| + 2k = m' . By removing edges, 
we would construct a proper subgraph of h with the same properties 
as h. 

To build g, we complete h in both cases by the part of g not concerned 
by the renormalization. □ 

Lemma 3.4. Let g = (V g ,E g ) be a (k, I) -redundant subgraph of G, 
minimal (in the sense that no proper subgraph of g is (k, I) -redundant). 
Then i) or ii) is true: 

i) g is a subset of H (that is g is included in the part of G to be 
renormalized). 

ii) g = R{g) is a (k, I) -redundant subgraph of G = R(G) (that is the 
image of g by renormalization is still (k, I) -redundant). 

Proof: Suppose i) is not true, and A,B e g. Let us call h = (Vh, Eh) 
the intersection of g and H, and h = R(h). We call fh the number of 
edges of h. The set of vertices and edges of g are called V g and E g 
respectively. Then \E g \ = \E g \ — fh + \E h \ and \V g \ = \V g \ + \Vh\ — 2. 
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Thus 

(3.5) \E g \ - k\V g \ + I = \E~ g \ - k | Vg | + I + \E h \ - k\V h \ + 2 jfe 

Now, /i is (k,l)-sparse since g is a minimal redundant graph. Thus 
h G Sab (fi) j and 

(3.6) m> - jfe|Vft| + 2A; 



Putting together Eqs. (3.5) and (3.6), an using \E g \ — k\V g \ + 1 > we 
get 

(3.7) \E~ g \ - k\V s \ +1 > 

This proves that ii) is true. 

If i) is not true and neither A nor B are in g, then ft, is actually 
empty: otherwise, g would be a minimal redundant subgraph with two 
disconnected components, which is impossible. Then the renormaliza- 
tion does not modify g and ii) is true. 

If i) is not true and A or B, but not both, are in g, then either h has 
no edge, and we are done since g is not modified by renormalization, 
or g is a minimal redundant graph made of two components that share 
only one vertex. It is easy to show by enumerating vertices and edges 
that this is impossible as soon as I > k. This ends the proof. □ 

NB: Notice that for < I < k, this lemma is not true. Fig. [5] shows 
a counter example for (k,l) = (2,1). It is not clear to us how it is 
possible to define a useful renormalization transformation in this case. 

Lemma 3.5. Let g = (v,e) be a (k, I) -redundant subgraph of G. Then 
there exists g, (k, I) -redundant subgraph of G such that g = R(g) (that 
is: g is the image by renormalization of a (k, I) -redundant subgraph of 
G). 

Proof: If e does not contain any A — B edge, then it is enough to 
take g = (v, e). 

We assume now that e contain m > edges linking A and B. This 
implies that Mab(H) > (since Mab(H) is the number of A — B edges 
in G). 



Using the same reasoning as in Lemma 3.3 , we construct h = (Vh,Eh) 
an element of Sab{H) such that \Eh\ — k\Vh\ + 2k = rh and such that 
the renormalization of h produces exactly fh edges between A and B. 

Define now g = (V g , E g ), with V g — v U Vh and E g contains E h and 
e minus the fh edges between A and B. Then R(g) = g, and 

\E g \-k\Vg\+l = \e\—rh+\E h \-k\v\+2k-k\Vh\+l = \e\-k\v\+l > 
Then g is redundant. □ 
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Figure 5. For (k,l) = (2, 1), the subgraph H (induced 
by vertices (A, B, X, Y)) does not contain any redundant 
edge, and is renormalized into a subgraph with a single 
edge between A and B. The renormalized subgraph is 
sparse. However, the subgraph induced by the vertices 
(A, X, Y, Z, T) is redundant. 



Proof of |3.2| i): 

Rigid connection in G =>- rigid connection in G: Suppose A and B are 
(k, I) -rigidly connected in G. Then there exists G\ = (Vi,Ei) a (k, l)- 
tight subgraph of G containing A and B. Let us call H\ = (V^jEh-J 
the intersection of G\ and the subgraph to be renormalized H; G± = 
R(Gi) = (Vi, Ei), Hi = R(Hi) = (V Hl ,E Hl ) the subgraphs of G images 
of Gi,Hi by renormalization; Hi has only two vertices A and B, and 
a certain number of edges between them. 
We have 

|Vi| = |V X | + \V Hl \-2- | | = |^i| + \E Hl \ - \E Hl \ 

Hi is a subgraph of Gi which is tight, so that Hi is sparse. The 
subgraph of Gi defined by the vertices Vi \V/y 1 U {A, B} is also sparse. 
This yields the two inequalities 

k\V Hl \ -2k< \E Hl \ < k\V Hl \ - I 

Then, by renormalization of Hi, we have 

1^1 = 1^1-^1^1+2* 

From this and the fact that Gi is tight, we obtain finally \ Vi \ = k\Ei\ —I. 

Now we want to prove that Gi is sparse. Let g = (v, e) be a s ubgr aph 
of Gi. We have to show that \v\ — k\e\ + I < 0. By Lemma 



3.3 



we 



construct g = (v , e). h — (Vh, E h ) is the intersection of g with Hi, and 



14 



J. BARRE 



m is the number of edges in the renormalization of h. Then using point 



Hi) of lemma 3.3 and 



\v\ = \v\ + \Vh\ — 2 ; |e| = |e| — in + \E h \ , 

we conclude 

|e| - k\v\ +1 < 

This proves that G\ is (k, Z)-tight, so that A and B are rigidly con- 
nected in G. 

Rigid connection in G =>- rigid connection in G: We assume now 
that A and B are rigidly connected in G, and want to show this is true 
also in G. Let G\ = (Vi, Ei) be a (k, Z)-tight subgraph of G containing 



A and B. Using Lemma 3.3, we construct G\ = (Vi,Ei) subgraph 
of G, and call Hi = (Vh^E^) the intersection of G\ with H. First, 
reasoning as above, we have \E\\ = A;|Vi| — I. 

We now have to prove that G\ does not contain any (k, /)-redundant 
subgraph. Suppose it is not the case, and call g\ a minimal redun- 



dant subgraph of G\. Then we have the alternative of Lemma 3.4| i) 



cannot be true because g\ would be a subgraph of Hi, which is sparse 



by Lemma 3.3 ii) cannot be true either, because we would have con- 
structed a redundant subgraph of Gi, which is tight by hypothesis. We 
conclude that Gi is tight, and that A and B are rigidly connected in G\. 



Proof of 3.2 ii): 



We want to prove the formula 

(3.8) n red (G) = n red (G) + n red (H) 

To prove this, we follow the algorithm defining n red (G), as described in 
Sec. [2] Starting from the vertices V and no edge, we add the edges in E 
one by one, starting with the edges which are in H, the subgraph to be 
renormalized. If the newly added edge is redundant, we discard it and 
add one to the count of n red {G). This way, we construct a sequence of 
sparse graphs 

Go = (V, 0), Gi, . . . , Gm h , • • • , Gm = (V, Em) , 

where M H = \E H \ and M = \E\. We have n red (G) = M - \E M \. E M 
contains the edges in E which have not been found redundant in the 
sequential edge addition process. 

Consider the associated sequence of renormalized graphs 

Gq = R(Gq), Gi, . . . , Gm h , ■ ■ ■ , Gm = R(Gm) 
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This is a sequence of sparse graphs thanks to lemma 3J3 and along 
this sequence, all edges in E are added one by one. We will show we 
can use this sequence to count the number of redundant edges in G. 

Adding the edges of H: 

When adding one by one the edges of H, exactly n re d(H) edges are 
discarded, and are added to the count of n re d{G). After addition of 
all edges of H, Gm h contains exactly Mab(H) edges between vertices 
A and B. Along this sequence of edges additions, (G n ) n <M H contains 
only a number smaller than 2k — I of edges connecting A and B, and 
no other edge. Thus, the count of n re d{G) remains 0. 

Adding an edge e E E\Eh (ie an edge in G, but not in H) : 

We start with the graphs G n = (V,E n ) and G n = (V,E n ), with 
n > Mjj. In this case, trying to add an edge in G corresponds to 
trying to add the same edge in G. We will show that e makes either 
both G n and G n redundant, or none of the two. 

Case 1 : (V, E n U {e}) is not redundant: the added edge is accepted, 
G n+ i = (V, E n U {e}), and the count of n re d{G) is not modified. 
Then (V,E n U {e}) is not redundant: otherwise, by lemma 3.5, it 
would be the image by renormalisation of a redundant subgraph of 
(V, E n U {e}), which is impossible. Then G n+ i = (V,E n U {e})), and 
neither the count of n re d{G) nor the count of n re d{G) are modified. 

Case 2: (V,E n U {e}) is redundant: the added edge is discarded, 
G n +i = (V, E n ) = G n , and the count of n re d{G) is increased by 1. 
Then (V, E n U{e}) is redundant. Indeed, let g by a minimal redundant 
subgraph of (V,E n U {e}). It cannot be included in H, as it contains 



the edge e. Then, by lemma 3.4, R(g) is a redundant subgraph of 



(V,E n U {e}). Thus the count of n re d(G) is also increased by 1, and 
Repeating this until all edges in E\E^ have been added proves for- 



mula (3.8). □ 

Notice that formula (3.8) is wrong for the graph of Fig. [5j with 
(k,l)-(2,\). 



3.3. Hierarchical graphs. If it is possible to iterate the renormal- 
ization transformation we have just defined until the graph becomes 
trivial, then the problems of percolation and counting the number of 
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Figure 6. A hierarchical graph constructed from a 
"Wheatstone bridge" elementary cell. From left to right, 
the level is t — 0, t — 1, t = 2. This is the graph which 
is used in Section |4j where in addition each edge may 
have multiplicity up to 3. 




redundant edges are exactly solved. It is indeed possible to define 
some graphs on which this procedure can be carried out completely, 
thus providing exactly solvable models of (k, I) percolation that go be- 
yond trees and random graphs. These graphs are called "hierarchical 
graphs" [19J, and are defined as follows. We start from two vertices, 
connected by one edge. The graph is then constructed iteratively; at 
each step, all edges are replaced by a given elementary cell. From each 
type of elementary cell, one thus constructs a hierarchical graph. A 
graph where this replacement procedure has been iterated t times will 
be called a level t hierarchical graph. An example is given on Fig. [6j 

Consider now within a hierarchical graph one elementary cell be- 
tween vertices A and B. It clearly satisfies the renormalizability condi- 



tion of Sec ]3.2[ all edges going outside of the cell are connected to the 
outer vertices A and B, whereas the inner vertices are only linked within 
the cell. Thus, the renormalization procedure described in Sec. 3.2 is 
exact when applied to an elementary cell. Furthermore, if the renor- 
malization procedure is applied to all elementary cells of a level t hier- 
archical graph, the resulting renormalized graph is again a hierarchical 
graph, of level t — 1. These remarks explain why the renormalization 
procedure allows to solve exactly (k, /)-rigidity models on hierarchical 
graphs. 

4. An exactly solved model 

In this section, we apply the formalism developed in the previous 
section to a body-bar rigidity model: that is we study (3, 3) percolation 
on a multigraph G, with N vertices and M edges. The multigraph G is 
constructed starting from the "Wheatstone bridge" hierarchical lattice, 
described on Fig. |6j 

We define the percolation problem as follows: each edge in the hi- 
erarchical lattice has multiplicity 0, 1, 2 or 3 with probability Po,Pi,P2 
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and pz = 1 — po — p± — p 2 (an edge with multiplicity is absent). 
Notice that if edges are either absent or have multiplicity 3 (that is 
Pi = Pi = 0), a subgraph is (3,3) rigid if and only if it is connected. 
In other words, (3,3) rigidity with triple edges is equivalent to (1,1) 
rigidity, that is ordinary percolation, with simple edges. For the special 
values of the parameters p\ = p 2 = 0, this model thus contains ordinary 
percolation; since (3, 3) rigidity percolation is supposed to belong to a 
different universality class, it is interesting to study the whole phase 
diagram of this model. 

Let us start from a level t hierarchical lattice, with t large. The 
renormalization rule of Sec. [3] applied to each unit cell transforms the 
level t lattice to a level t — 1 lattice, and induces a transformation of 
the parameters pf. 

(Pa ,P* ,P2 ,P3 ) = <p[(Pa\Pi \pi jPs)] 

From an analysis of the renormalization transformation, it is not diffi- 
cult to obtain the explicit expression of <p, but it is very tedious. The 
details are given in the Appendix. Assuming that we start from a large 
graph with t ^> 1, we iterate this renormalization transformation. To 
understand the phase diagram of the model, one needs now to study 
the renormalization flow induced by tp; of particular interest are the 
fixed points of (p. Note that the four dimensional space of parameters 
is actually easily reduced to three dimensions, since P0+P1+P2+P3 = 1- 
By inspection of the complicated expression for <p, three fixed points 
are easily found: 

l ) P*flo P = (Po = 1, Pi = 0, p 2 = 0, p 3 = 0). This is the empty graph 
fixed point, corresponding to the floppy phase. 

w) P*ri g = (Po = 0, pi = 0,p 2 = 0,p 3 = 1). This is the full graph 
fixed point, corresponding to the rigid phase. 

Hi) p* perc = (po = 0.5, pi = 0,p 2 = 0,p 3 = 0.5). This is a criti- 
cal fixed point, corresponding to ordinary percolation (because in this 
case, edges are either absent, or have multiplicity 3). 

Using a Newton- Raphson scheme and scanning the whole parameter 
space, we have found another fixed point: 

w ) P*crit = (Po ^ 0.326, pi ~ 0.174, p 2 ~ 0.174, p 3 ~ 0.326). This 
fixed point may be found also by noting that the surface pi — p 2 , po = 
P3 is stable by ip. Using the normalization condition P0+P1+P2+P3 = 1, 
looking for a fixed point on this surface is then a one dimensional 
problem. We have not found any other fixed point in the domain 

Pi >0,P0+Pl+P2+P3 = 1- 
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The trivial fixed points p* flop and p* rig are stable. The fixed point 
P* P erc nas three unstable directions. The fixed point p* crit has only 
one unstable direction. The renormalization flow is then as follows: 

• The three dimensional parameter space is divided by a criti- 
cal hypersurface containing p* perc and p* crit . On one side of 
the surface, the renormalization flows approaches the empty 
("floppy", non percolating) fixed point; on the other side, it 
approaches the full fixed point ("rigid", percolating). 

• On the critical hypersurface, the flow is attracted by p* crit , the 
"rigidity percolation" critical fixed point. 

We conclude that except for very special choices of parameters p\ = 
P2 = 0, the large scale critical properties of this model are described 
by the "rigidity percolation" critical fixed point. 

5. Conclusion 

We have introduced and rigorously justified a renormalization trans- 
formation adapted to the study of (k, /)-rigidity, for k < I < 2k, which 
generalizes the well known procedure for ordinary percolation, and the 
procedure for bar-joint rigidity used in [16]. This method allows to 
solve exactly (k, /)-percolation problems on hierarchical graphs. We 
have provided such an example, which has the interesting feature of 
showing both "ordinary percolation" and "rigidity percolation" behav- 
ior for different values of the parameters. 

Rather than solving exactly problems on hierarchical graphs, such 
renormalization transformations might be used to provide approxi- 
mate solutions for problems on more realistic 2D or 3D lattices. This 
work would then provide an approximate analytical tool in the gen- 
eral study of (k, Z)-percolation problems and their universality classes, 
a field which is still widely open. 

This work is supported by the ANR-09-JCJC-009401 INTERLOP 
project. 

6. Appendix 

We give here the details of the computations yielding the renormal- 



ization function ip. The problem is simple: apply the rule of Sec. 34^ to 
an elementary "Wheatstone bridge" cell. Since there are 5 edges, and 
each edge may be absent, single, double or triple, there are 4 5 = 1024 
different configurations of which we have to study the renormalization 



according to the rule of Sec. 3.1 The large number of configurations 



to enumerate is the only difficulty Prop. 3.1 is useful to simplify these 



1!) 



computations, as it allows to remove all redundant edges from the be- 
ginning. 

In the following tables, we classify the configurations according to 
the multiplicity of their edges. For instance, an edge multiplicity 32210 
means that one edge has multiplicity 3, two have multiplicity 2, one has 
multiplicity 1 and one is absent. Clearly, depending on how these edges 
are distributed on the "Wheatstone bridge", the number of renormal- 
ized edges may be different. The second column of the tables contains 
the number of edges of the renormalized cell, and the third column 
is the combinatorial factor corresponding to the number of configu- 
rations with the given edge multiplicity yielding the given number of 
renormalized edges. 

We have collected the results in 6 tables, according to the number 
of edges with multiplicity 3. 



5 edges with multiplicity 3 



edges multiplicity in H 


nbr of renormalized edges 


nbr configurations 


33333 


3 


1 


4 edj 


^es with multiplicity 3 


edges multiplicity in H 


nbr of renormalized edges 


nbr configurations 


33332 
33331 
33330 


3 
3 
3 


5 
5 
5 



3 edges with multiplicity 3 



edges multiplicity in H 


nbr of renormalized edges 


nbr configurations 


33322 


3 


10 


33321 


3 


20 


33320 


3 


16 


33320 


2 


4 


33311 


3 


8 


33311 


2 


2 


33310 


3 


16 


33310 


1 


4 


33300 


3 


8 


33300 





2 



20 
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2 edj 


ges with multiplicity 3 




1 1 J_ * 1 * * J_ * TT 

edges multiplicity m H 


1 £ 1*11 

nbr of renormalized edges 


nbr configurations 


33222 


3 


10 


33221 


3 


30 


33220 


3 


18 


33220 


2 


12 


33211 


3 


24 


33211 


2 


6 


33210 


3 


36 


33210 


2 


12 


33210 


1 


12 


33200 


3 


6 


33200 


2 


18 


33200 





6 


33111 


3 


6 


33111 


2 


4 


33110 


3 


6 


qqi 1 n 

OOl 1U 


o 
A 


1 


33110 


1 


12 


33100 


3 


6 


33100 


1 


18 


33100 





6 


33000 


3 


2 


33000 
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Collecting the information from these tables, we obtain the expres- 
sion for the renormalization function (p. Calling the renormalized prob- 
abilities p' Q , p[, p 2 , p 3 , we have 

p' 3 = =pl + 5pjp + 8plp 2 + 2p 2 p 3 +p 5 2 + 5p 4 2Pl + 6p 2 3 p lP 2 + 6p 2 3 p 2 lPo + ±p\p\ 
+ lQplpiPo + 8p 3 3 pj + 5pgj9i + 5p 3 p 2 + 8p 3 p 3 2 p + 6p 2 3 p 2 p 2 + ^pIpIpq 
+lpplpi + IOP3P2 + I6P3P2P0 + 5p^ 2 + 20p^ 2 p! + 3()p£p!pi + 24 P3 Wi 
+36j93j9 2 piPo + 20^3^1 + 21p$p\p\ 

p' 2 = lOplpj + 12plplp + 6p|p? + 2p\p\ + 4p 3 p 2 Po + ^PspIpI + ^P3pIpo 
+ I8p 2 3 p 2 p 2 + ^plp2Po + §p\pip\ + I2p 2 3 p 2 pip + Qp 3 p 2 2 p\ + <±8pzp 2 2 pipo 
+20p 3 p 2 pl + 12p 3 p 2 plpo + I2p 3 p 2 pip 2 + 5p 2 p + 12p 2 3 p 2 2 p 

Pi = ^pIpI + QpIpI + ^pIpI + QpIpIpo + Qp 2 2 PiP 2 o + 2OP2P1P0 + 4j9 3 PiPo + 12 P3P 2 iP 2 o 
+12p 3 pip + 5p 3 pi + 18p 3 piPo + 12 PlPiPo + ^pIpiPo + ^Pzp\p\ + V2p 2 3 p 2 p x p$ 
+ I2p 3 p 2 2 pip + A8p 3 p 2 p\p + \2p 3 p 2 p x p\ 

Po = Pi + SpsPo + 8P3P0 + 2p 3 Po + 5p 2 Po + %pIpI + 4 ^2Po + Pi + 5PiPo + 10p?Po 
+ 10p^o + 5 PiPo + 5 P2Pi + 20p 2 plp + 30p 2 plpl + 20p 2 p 1 pf ) + 2Ap 2 2 p\po 
+2Ap 2 2 pip 2 + lQp 3 pipl + I8p 3 plp 2 + 8p 3 p\p Q + Qp 2 3 pipl + 36p 3 p 2 PiPo 
+6P3P2P0 + §p\pip\ + 16p 3 p 2 Po 
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1 ed, 



ge with multiplicity 3 



edges multiplicity in H 


nbr of renormalized edges 


nbr configurations 


32222 


3 


5 


32221 


3 


20 


32220 


3 


8 


32220 


2 


12 


32211 


3 


24 


32211 


2 


6 


32210 


2 


48 


32210 


1 


12 


32200 


2 


12 


32200 


1 


12 


32200 





6 


32111 


2 


20 


32110 


2 


12 


32110 


1 


48 


32100 


2 


12 


32100 


1 


12 


32100 





36 


32000 


2 


4 


32000 





16 


31111 


1 


5 


31110 


1 


12 


31110 





8 


31100 


1 


12 


31100 





18 


31000 


1 


4 


31000 





16 


30000 
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No edge with multiplicity 3 



edges multiplicity in H 


nbr of renormalized edges 


nbr configurations 


22222 


3 


1 


22221 


3 


5 


22220 


2 


5 


22211 


2 


10 


22210 


1 


20 
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1 


6 
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